Search for: All records

Lyman α emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to overdense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assumed to mark overdense, ionized bubbles. Here, we develop the first framework to quantify the ionization and density fields of high-z galaxy associations. We explore the interplay between (i) the large-scale density of a survey field, (ii) Poisson noise due to the small number density of bright sources at high redshifts (z ∼ 7), and (iii) the effects of the ionized fraction on the observation of LAEs. We use Bayesian statistics, a simple model of reionization, and a Monte Carlo simulation to construct a more comprehensive method for calculating the large-scale density of LAE regions than previous works. We find that Poisson noise has a strong effect on the inferred density of a region and show how the ionized fraction can be inferred. We then apply our framework to the strongest association yet identified: Hu et al. found 14 LAEs in a volume of ∼50 000 cMpc3 inside the COSMOS field at z ∼ 7. We show that this is most likely a 2.4σ overdensity inside of an ionized or nearly ionized bubble. We also show that this LAE association implies that the global ionized fraction is $\bar{Q} = 0.59^{+0.10}_{-0.11}$, within the context of a simple reionization model.

 

Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization. As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch and could thus potentially pinpoint some of the earliest protoclusters. In this Letter, we report the discovery of a massive protocluster around a luminous quasar atz= 6.63. This protocluster is anchored by the quasar and includes three [Cii] emitters atz∼ 6.63, 12 spectroscopically confirmed Lyαemitters (LAEs) at 6.54 <z≤ 6.64, and a large number of narrow-band-imaging selected LAE candidates at the same redshift. This structure has an overall overdensity of$\delta ={3.3}_{-0.9}^{+1.1}$within ∼35 × 74 cMpc²on the sky and an extreme overdensity ofδ> 30 in its central region (i.e.,R≲ 2 cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of${6.9}_{-1.4}^{+1.2}\times {10}^{15}{M}_{\odot }$at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE, which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium.

 

Observed scatter in the Lyαopacity of quasar sightlines atz< 6 has motivated measurements of the correlation between Lyαopacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Lyαtroughs atz∼ 5.7 with a deficit of Lyαemitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift toward thez= 6.02 quasar SDSS J1306+0356 and thez= 6.17 quasar PSO J359-06. We find that both fields are underdense in LAEs within 10h⁻¹Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Lyαtroughs. We combine our observations with three additional fields from the literature and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in intergalactic medium (IGM) Lyαopacities. There is tension at low opacities, however, as the models tend to associate lower IGM Lyαopacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent.

 

Fluctuations in Lyman-α (Ly α) forest transmission towards high-z quasars are partially sourced from spatial fluctuations in the ultraviolet background, the level of which are set by the mean free path of ionizing photons (λmfp). The autocorrelation function of Ly α forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to λmfp. Recent measurements at z ∼ 6 suggest a rapid evolution of λmfp at z > 5.0 which would leave a signature in the evolution of the autocorrelation function. For this forecast, we model mock Ly α forest data with properties similar to the XQR-30 extended data set at 5.4 ≤ z ≤ 6.0. At each z, we investigate 100 mock data sets and an ideal case where mock data matches model values of the autocorrelation function. For ideal data with λmfp = 9.0 cMpc at z = 6.0, we recover $\lambda _{\text{mfp}}=12^{+6}_{-3}$ cMpc. This precision is comparable to direct measurements of λmfp from the stacking of quasar spectra beyond the Lyman limit. Hypothetical high-resolution data leads to a $\sim 40~{{\ \rm per\ cent}}$ reduction in the error bars over all z. The distribution of mock values of the autocorrelation function in this work is highly non-Gaussian for high-z, which should caution work with other statistics of the high-z Ly α forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements.

 

Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21-cm power spectrum. One of these parameters, Rmax, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum distance that ionizing photons can travel. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically motivated prescriptions of the mean free path (MFP) of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in seminumerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly defined MFP prescriptions more effectively suppress large-scale structures in the ionization field in seminumerical reionization simulations compared to the standard Rmax approach, and the magnitude of the MFP modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcoming 21-cm power spectrum observations. Finally, we show that in our model, the MFP plays a larger role in regulating the reionization history than in models using Rmax.

 

The mean free path of ionizing photons,λ_mfp, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements ofλ_mfpfrom QSO spectra over the redshift range 5 <z< 6, including the first measurements atz≃ 5.3 and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs nearz∼ 5.5, for which we also acquire new [Cii] 158μm redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find${\lambda }_{\mathrm{mfp}}={9.33}_{-1.80}^{+2.06}$,${5.40}_{-1.40}^{+1.47}$,${3.31}_{-1.34}^{+2.74}$, and${0.81}_{-0.48}^{+0.73}$pMpc atz= 5.08, 5.31, 5.65, and 5.93, respectively. Our results demonstrate thatλ_mfpincreases steadily and rapidly with time over 5 <z< 6. Notably, we find thatλ_mfpdeviates significantly from predictions based on a fully ionized and relaxed IGM as late asz= 5.3. By comparing our results to model predictions and indirectλ_mfpconstraints based on IGM Lyαopacity, we find that the evolution ofλ_mfpis consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift 6.

 

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low (${\overline{x}}_{\mathrm{H}\mathrm{I}}\lesssim {10}^{-3}$) or close to unity (${\overline{x}}_{\mathrm{H}\mathrm{I}}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on${\overline{x}}_{\mathrm{H}\mathrm{I}}$atz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.3)<0.79\pm 0.04$(1σ),${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.5)<0.87\pm 0.03$(1σ), and${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.7)<{0.94}_{-0.09}^{+0.06}$(1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

 

We present measurements of black hole masses and Eddington ratios (λ_Edd) for a sample of 38 bright (M₁₄₅₀< −24.4 mag) quasars at 5.8 ≲z≲ 7.5, derived from Very Large Telescope/X–shooter near–IR spectroscopy of their broad Civand Mgiiemission lines. The black hole masses (on average,M_BH∼ 4.6 × 10⁹M_⊙) and accretion rates (0.1 ≲λ_Edd≲ 1.0) are broadly consistent with that of similarly luminous 0.3 ≲z≲ 2.3 quasars, but there is evidence for a mild increase in the Eddington ratio abovez≳ 6. Combined with deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 158μm line from the host galaxies and VLT/MUSE investigations of the extended Lyαhalos, this study provides fundamental clues to models of the formation and growth of the first massive galaxies and black holes. Compared to local scaling relations,z≳ 5.7 black holes appear to be over-massive relative to their hosts, with accretion properties that do not change with host galaxy morphologies. Assuming that the kinematics of theT∼ 10⁴K gas, traced by the extended Lyαhalos, are dominated by the gravitational potential of the dark matter halo, we observe a similar relation between black hole mass and circular velocity as reported forz∼ 0 galaxies. These results paint a picture where the first supermassive black holes reside in massive halos atz≳ 6 and lead the first stages of galaxy formation by rapidly growing in mass with a duty cycle of order unity. The duty cycle needs to drastically drop toward lower redshifts, while the host galaxies continue forming stars at a rate of hundreds of solar masses per year, sustained by the large reservoirs of cool gas surrounding them.

 

Abstract The observed large-scale scatter in Ly α opacity of the intergalactic medium at z < 6 implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly α opacity and density. We present selections of z = 5.7 Ly α -emitting galaxies (LAEs) in the fields surrounding two highly opaque quasar sightlines with long Ly α troughs. The fields lie toward the z = 6.0 quasar ULAS J0148+0600, for which we reanalyze previously published results using improved photometric selection, and toward the z = 6.15 quasar SDSS J1250+3130, for which results are presented here for the first time. In both fields, we report a deficit of LAEs within 20 h −1 Mpc of the quasar. The association of highly opaque sightlines with galaxy underdensities in these two fields is consistent with models in which the scatter in Ly α opacity is driven by large-scale fluctuations in the ionizing UV background or by an ultra-late reionization that has not yet concluded at z = 5.7.